Effect of Noise on DNA Sequencing via Transverse Electronic Transport 



Matt Krems 1 , Michael Zwolak 2 , Yuriy V. Pershin 3 , Massimiliano Di Ventra 1 

1 Department of Physics, University of California, San Diego, La Jolla, CA 92093 
2 Theoretical Division, MS-B213, Los Alamos National Laboratory, Los Alamos, NM 87545 
3 Department of Physics and Astronomy and USC Nanocenter, 
University of South Carolina, Columbia, SC 29208 
(Dated: June 9, 2009) 

Previous theoretical studies have shown that measuring the transverse current across DNA strands 
while they translocate through a nanopore or channel may provide a statistically distinguishable 
signature of the DNA bases, and may thus allow for rapid DNA sequencing. However, fluctuations 
of the environment, such as ionic and DNA motion, introduce important scattering processes that 
may affect the viability of this approach to sequencing. To understand this issue, we have analyzed 
a simple model that captures the role of this complex environment in electronic dephasing and its 
ability to remove charge carriers from current-carrying states. We find that these effects do not 
strongly influence the current distributions due to the off-resonant nature of tunneling through the 
nucleotides - a result we expect to be a common feature of transport in molecular junctions. In 
particular, only large scattering strengths, as compared to the energetic gap between the molecular 
states and the Fermi level, significantly alter the form of the current distributions. Since this gap 
itself is quite large, the current distributions remain protected from this type of noise, further sup- 
porting the possibility of using transverse electronic transport measurements for DNA sequencing. 



Introduction 

The prospect of sequencing an entire human genome 
for less than $1000 USD in a matter of hours is 
becoming closer to reality [TJ [U [3J. The original 
DNA-nanopore experiments of Kasianowicz et al. jl] 
showed polynucleotides can be pulled through nanoscale 
pores and their translocation detected by measuring 
the consequent blockage of the ionic current through 
the pore. Since then, numerous experimental studies 
have been performed using biological [5j [5) [51 M 
HO] and synthetic nanopores QTJ E2 [T3J El E33 US] 
which probe various physical properties of translocating 
polynucleotides. This has fueled an enormous amount 
of research into novel sequencing proposals based on 
nanopores or nanochannels [TJ [H [3J . 

One sequencing idea suggests detecting transverse elec- 
tron currents as single-stranded DNA (ss-DNA) translo- 
cates through a pore [T5J EE HTJ1 121] ■ Previous theo- 
retical work showed the four DNA nucleotides possess 
statistically distinguishable electronic signatures [TH [T3] 
in the form of current distributions when accounting for 
structural distortions and partial control of the DNA dy- 
namics (i.e., by a transverse field) [THJ [201 [21]. These 
results indicate DNA sequencing is, in principle, possi- 
ble via transverse current measurements. However, such 
studies have neglected scattering processes, such as fluc- 
tuations of the environment, which introduce electronic 
noise, and may thus affect the ability to distinguish the 
bases. 

Recently, experimentalists have successfully embed- 
ded electrodes into solid state nanopores and nanochan- 
nels (221 [23j [211 [25] and are getting closer to measuring 
electronic currents with single nucleotides present in the 




FIG. 1: Schematic representation of ss-DNA translocating 
through a pore while the transverse electronic current is col- 
lected. The light (purple) atoms are the silicon nitride pore 
and the dark (black) atoms represent the electrode surfaces 
within our molecular dynamics simulations. The single strand 
of DNA translocates through the pore pulled by a longitudi- 
nal electric field, E\\ , and the nucleotides also experience a 
transverse electric field, Ej_. The white arrows around the 
DNA base indicate an acoustic phonon-like motion which 
contributes to the noise. The visualization was made with 
VMD [T7|. 



gap between the electrodes. When the latter is achieved, 
one question which will arise is how does the noise in- 
duced by the environment - noise which is beyond that 
due to "static" structural distortions of the nucleotides - 
affects the nucleotides' electronic signatures, i.e., the cur- 
rent distributions. The environment is composed of ionic 
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and water fluctuations and other excitations which may 
drastically affect the electron dynamics, and thus the cur- 
rent |26j . To complicate matters, the liquid environment 
can scatter electrons out of their current-carrying states 
by absorbing them into the solution and allowing the lon- 
gitudinal field (that pulls the DNA through the pore) to 
carry them away. The influence of these and related fac- 
tors can be very important, as seen in previous studies 
of electronic transport through DNA [S7J HHJ , and so 
far no study has examined such effects in detail. 

In this article, we address these issues theoretically. 
Clearly, a fully time-dependent calculation with inclu- 
sion of all these types of scattering processes would be 
ideal [26]. However, the complexity of the problem we 
consider, both in the number of atoms involved and the 
type of scattering processes to take into account, makes 
this type of dynamical calculation unrealistic at present. 
Instead, we use a simplified model to capture some of the 
physics we deem important and leave a time-dependent 
treatment for future investigation. 

In general, one expects any type of electronic noise 
to eventually destroy the capability to distinguish the 
DNA bases once its strength is sufficiently large. Indeed, 
we do find this type of behavior. However, the noise 
strength at which the electronic transport is negatively 
influenced is very large, beyond the strength one would 
expect in realistic experimental situations. This is due 
to the off-resonant nature of tunneling through the nu- 
cleotides, and we thus expect this result to be a common 
feature of molecular junctions. In other words, the sep- 
aration of the energy levels of the nucleotides from the 
equilibrium Fermi level "protects" the electronic signa- 
ture of the bases. The present study will thus help re- 
searchers understand future experimental data, and pro- 
vides further support to the viability of DNA sequencing 
via transverse electronic transport. 

Setup and Methods 

As our starting point, we employ molecular dynamics 
simulations performed with NAMD2 [30J to pull homoge- 
neous ss-DNA through a SiaN4 nanopore with embedded 
gold electrodes. Our basic setup is shown in Fig. [T] and 
is the same as that used in previous work [THl [2U], except 
the new trajectories here correspond to longer simulation 
times. These trajectories give us the real-time atomistic 
structure of ss-DNA as it propagates through the pore. 
With these structures, we calculate the electronic trans- 
port in the transverse direction across the pore. In the 
latter calculations, we include the effect of noise as dis- 
cussed below. 

The details of the simulations are as follows. The pore 
is made of 2.4 nm thick silicon nitride material in the 
/3-phase. The nanopore hole has a double conical shape 
with a minimum diameter of 1.4 nm located at the center 



of the membrane and an outer diameter of 2.5 nm (see 
Fig. [TJ. The inner diameter is chosen wide enough such 
that ss-DNA is able to pass through but narrow enough 
that an appreciable tunneling current can be detected. 
The nanopore is then solvated in a TIP3 water sphere 
of 6.0 nm radius with spherical boundary conditions in 
an NVT ensemble and with a 1 M solution of potassium 
and chlorine ions. The CHARMM27 force field [HI [32] is 
used for the interaction of DNA, water, and ions, while 
UFF [33] parameters are used for the interaction of the 
SiaN4 membrane and other atoms. The SiaN4 atoms are 
assumed to be fixed during the simulation (this does not 
affect the conclusions). A 1 fs timestep is used and the 
system temperature is kept at room temperature with a 
Langevin dampening parameter of 0.2 ps -1 in the equa- 
tions of motion [34 . The van der Waals interactions are 
gradually cut off starting at 10 A from the atom until 
reaching zero interaction 12 A away. The energy was 
initially minimized in 1000 time steps. 

A single strand of DNA is constructed by removing 
one strand from a helical, double-stranded polynucleotide 
created using the Nucleic Acid Builder of the Amber Tools 
package [35] • At the initial time of the simulation, the 
ss-DNA is placed parallel to the pore axis with the first 
base just inside the pore. The ss-DNA is driven through 
the pore with a global electric field of 6 kcal/ (mol Ae) to 
achieve reasonable simulation times. In the calculation 
of the electronic transport, the longitudinal pulling field 
is turned off and a transverse field (of the same mag- 
nitude as that driving the current) is turned on at a 
moment when a base is between the electrodes. This 
approximates the situation when the transverse field is 
much larger than the longitudinal field. We envision this 
as the typical operating regime for a sequencing device 
as it allows for the suppression of a significant amount of 
structural distortion [19] . The particular time to stop the 
translocation is chosen by visual inspection. This stop- 
ping time is not particularly important because it gener- 
ally takes on the order of hundreds ps for the transverse 
field, E±, to align the nucleotide with the electrodes |20j . 
Single-stranded DNA differs from double-stranded DNA 
in that the persistence length of the polynucleotide is 
much shorter. This, in particular, allows for the base to 
quickly align with the perpendicular electric field. An 
example of this is reported in Fig. [2] where poly (C) 15 is 
such that a single C base is aligned parallel to a pair of 
opposite electrodes. A bias of 1 V oriented perpendicular 
to the base plane is then turned on. From the figure it 
is clear that, for this particular polynucleotide and ini- 
tial condition, the base and backbone reorient themselves 
towards the field within about 800 ps. This is also con- 
firmed by the currents as a function of time across two 
pairs of perpendicularly placed electrodes. At t — the 
largest current is from the pair of electrodes parallel to 
the plane of the base, while after 800 ps, the largest cur- 
rent is from the opposite pair of electrodes. It is also 
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evident from the figure that the rotation does not occur 
uniformly in time but it proceeds by fast rotations, fol- 
lowed by periods of time in which the system is temporar- 
ily trapped in a local energy minimum. Faster rotations 
have been observed with other initial conditions, trans- 
verse voltages and nucleotide strands [20] > but we cannot 
exclude the possibility that, for other initial conditions, 
longer times would be needed for a complete rotation of 
the bases. 
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FIG. 2: Currents as a function of time across two pairs of 
perpendicularly placed electrodes for poly (C) 15 with one base 
originally aligned parallel to a pair of opposite electrodes (see 
inset). The black current trace corresponds to the current 
from the black electrodes, and likewise the gray current trace 
corresponds to the gray electrodes. At time t = 0, a bias 
of 1 V oriented perpendicular to the base plane is switched 
on. The corresponding field aligns the base and backbone 
with the gray pair of electrodes (as shown in the inset), with 
a corresponding increase in the current across that pair of 
electrodes. 

The current calculations are performed within a single- 
particle scattering approach using a tight-binding Hamil- 
tonian (see, e.g., Ref. [26]). These calculations include 
water, although, within our approach, water has little di- 
rect effect on the current [20J ■ "Snapshots" of the atom- 
istic structure of ss-DNA between the gold electrodes are 
taken from the molecular dynamics at regular time inter- 
vals. These coordinate snapshots are used to obtain the 
tight-binding Hamiltonian. For each carbon, nitrogen, 
oxygen, and phosphorous atoms, s,p xi p y ,p z orbitals are 
used, whereas for gold and hydrogen only s orbitals are 
employed. The Fermi level is taken to be that of bulk 
gold 03]. 

To obtain the current across the ss-DNA, we use the 
retarded Green's function, 
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DNA 



(E) = 
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(1) 



ESdna — Hdna — S t — S h — E n 
where E is the energy, Sdna and Hdna are the overlap 



and Hamiltonian matrices, respectively, of the contents of 
the gap between the electrodes (we will call it electronic 
junction), St(t) are the self-energy terms associated with 
the interaction between the electrodes and the junction 
contents, and £„ is the self-energy associated with the 
noise. The Green's function for gold needed to calculate 
£ t (b) is approximated as in Ref. [3S]. We use a white- 
noise term, which corresponds to a noise timescale via 



n 



Im{£„}' 



(2) 



and we also take Re{£„} = (see discussion below). 
This timescale sets a decay time due to interaction with 
the environment. The latter can be thought of as a noise 
probe that interacts with the contents of the junction |261 

E2J. 

If we were to follow this type of reasoning we would 
then set the current in the probe to be zero and calculate 
the total transmission coefficient as 



T(E) = Tib (E) + T p (E) , 



where 



T tb (E) = Tr 



^tGDNA^bG 



t 

DNA 



(3) 



(4) 



is the transmission coefficient that directly couples elec- 
trodes that measure the current in the presence of the 
noise probe with r t (f,) = i (^t(b) ~ ^t(&))- This trans- 
mission contribution includes only elastic processes, as 
we discuss in more detail below. 
The other term is 



T P (E) 



TtnTnb 
Tn + T n b 



where 



Tfxu (E) — Tr r ^Gdna^vG\ 



DNA 



(5) 



(6) 



is instead the transmission from reservoir /i to v, namely 
it takes into account processes that can drive electrons 
out of the electrodes into the noise probe and vice versa. 
The current is then given by 



2e 



dET(E) [f t (E) - f b (E)} , 



(7) 



where f t n^ is the Fermi-Dirac function of the top (bot- 
tom) electrode [26j . The current distribution for a nu- 
cleotide is the distribution obtained from the various 
snapshots while the nucleotide fluctuates between the 
electrodes. 

We will later make a microscopic connection to the 
above transmission probability by starting with a Hamil- 
tonian for independent electrons coherently coupled to 
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FIG. 3: Transverse current versus time for poly(A)is at a 
transverse bias voltage of 1.0 V. The noise lowers the current 
slightly for r n = 10~ 15 s. Only at the unrealistic r n = 10 -16 s 
does the current shift significantly. Slower noise timescales 
give essentially the same current as the case with an infinite 
noise timescale. 



a phonon environment. However, this analysis leads us 
to conclude that in the complex liquid environment the 
term in Eq. [5] cannot correctly represent the physical sit- 
uation at hand. In fact, retention of such term would 
give rise to unrealistically large currents (several orders 
of magnitude larger than what we present here). While 
this result would naively suggest that such currents could 
in fact be observed in the present case, it is unlikely that 
the nanopore environment would allow for the coherent 
coupling between the electrons and excitations that gives 
this increased current. Furthermore, it is likely that due 
to the presence of the longitudinal field that drives the 
DNA through the pore the electrons scatter out of their 
current-carrying states. In this work we will then assume 
that current-carrying electrons can be scattered into the 
complex liquid environment and retain only the first term 
T t b in the transmission probability of Eq. [3] and analyze 
its effect as a function of the timescale strength r„ . This 
is equivalent to assuming that the liquid environment 
is represented by two probes connected to the junction, 
and the probes' electrochemical potentials are adjusted 
so that the combined current from the two probes into 
either electrode is zero. (45] 

Noise - As stated above, previous theoretical studies 
have shown the current distributions caused by DNA 
static structural distortions are statistically distinguish- 
able [H QH HH EU- These studies, however, have not 
included the effects of external noise. We focus specifi- 
cally on noise given by Eq. [2] because it represents many 
processes which happen in an experiment. These include 
fast processes, such as electronic interactions with bound 
waters or charges on the pore walls, and also slow pro- 



cesses, such as the dynamic movement of the DNA bases 
and ions. From visual inspection of the molecular dy- 
namics simulations, we observe that the bases fluctuate 
in a way reminiscent of acoustic phonons, i.e., we observe 
only low-energy excitations. An example of these exci- 
tations is represented in Fig. [I] where these slow oscilla- 
tions, while not periodic, are mostly in the longitudinal 
direction. No oscillations where the bases are, e.g., in 
a "breathing mode" , that is where the base itself is ex- 
panding and contracting, causing large energy relaxation, 
were observed. At low bias, these are also unlikely to be 
excited by the electrical current itself, so that we expect 
a low exchange of energy with the current-carrying elec- 
trons |38[ [39] . Furthermore, we assume the timescale 
for noise, Eq. [2j is a constant for all molecular states in 
the junction. In certain cases, this most likely overes- 
timates the effect of the noise, but, on the other hand, 
it misses "colored noise" effects, where, for instance, the 
noise has a strong component at a particular frequency. 
In the absence of a physical model for such noise which 
is supported by experiments, its effect is only specula- 
tive at this stage, and we thus defer its study for future 
research. 



Results and Discussion 

We have performed current calculations for 
some representative noise timescales [40 : r n = 
oo, 10~ 13 , 1CP 14 , 10 -15 , 10~ 16 s with transverse voltages 
of 0.1 V and 1.0 V. The timescale of 10~ 16 s is a particu- 
larly fast and unphysical timescale but was used to show 
the onset of major differences in the current and current 
distributions. 

For weak noise, (r„ = 10~ 13 s - 10 _14 s), the average 
current itself is essentially unchanged as well as the dis- 
tributions. The average percent change of an individ- 
ual current value for r n = 10~ 13 s is only about 0.1%. 
For r„ = 10~ 14 s, it is 1.5%. However, for a single 
current count, the current may vary by orders of mag- 
nitude due to the noise, further strengthening the ar- 
gument that a single base measurement is likely not 
enough to distinguish the bases [15]. From Figs. [3] and [4] 
T n = 10~ 15 s lowers the current on average and slightly 
alters the distributions. There is an average current re- 
duction of about 30%. At the unphysical fast timescale 
of 10~ 16 s, the current is significantly lowered and the 
distributions are pushed into an unmeasurable regime. 
However, we are not aware of a physical process that 
may cause such strong noise under the experimental con- 
ditions envisioned in this work. 

We have found above that even relatively strong noise 
does not negatively impact the current distributions. 
This may seem an unexpected result, and it will be help- 
ful for future experimental and theoretical efforts to un- 
derstand the reason for such an effect. We thus develop 
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FIG. 4: Probability distributions for poly(A)is with various 
noise timescales for a transverse bias voltage of 1.0 V. The 
very light dashed lines correspond to the bins used to pro- 
duce the current distributions. The solid lines are interpo- 
lated from the dashed ones. Like the current itself, only with 
very fast noise, T n — 10 _16 s, does the distribution change 
and shift appreciably. At t„ = 1CP 15 s, the distribution's 
mean shifts slightly and it broadens somewhat. 



a model system to understand this behavior, as well as 
the noise processes we are including. Our starting point 
is based on our previous work on transverse transport 
through DNA PO E3 EIfl E] • In an ideal configuration 
of a nucleotide between electrodes, the LUMO level of 
the base is closest to the gold Fermi level [TJ [18] and also 
couples well to both electrodes. Thus, it is reasonable to 
treat a nucleotide in the electronic junction as a single 
energy level, E . 

At this point we may consider a model Hamiltonian 
representing this level interacting with a bosonic envi- 
ronment 1461 



H = E dU+H de +H e +dUJ29k (b{ + b 



k 



(8) 

where cftd represents the occupation of the DNA LUMO 
level, i?de is the DNA-electrode interactions, and H e is 
the electrodes' Hamiltonian. The two remaining terms 
represent an interaction with a bosonic environment in 
the junction with interaction strength g k to each mode 
k. To get a tractable model, we make a few additional 
assumptions. First, we assume the junction DNA energy 
level is equally coupled to all levels of both electrodes and 
that we are at low enough bias and temperature (com- 
pared with electronic energies) that the electrodes band- 
width is effectively infinite. Second, we assume that the 
bosonic environment does not generate electronic corre- 
lations in the electrodes, which is reasonable for the small 
electrode coupling that we have here. Within these ap- 
proximations, the real-time retarded Green's function, 



Eq. [TJ becomes gTJ H2J 

Gdna (i) = -»0 (t) e-^'e^e-*® . (9) 

This Green's function includes the coupling to the elec- 
trodes through the factor e -7 *, where 7 is the coupling 
strength to both electrodes, and includes the coupling to 
the bosons through the factor e - ^'*) and the renormal- 
ized energy Eq. The bosonic term is 



tf> (t) = £ M- [n k (1 - e-** + K + 1) (1 - e— **)] , 
k ^ 

(10) 

where n k = 1/ (exp(/3o>fc) — 1) is the equilibrium occu- 
pation of mode k at inverse temperature (3. So long as 
the temperature is large compared to the boson cutoff 
frequency, uj c , then n k w l//3u! k and n k sw n k + 1, thus 
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In terms of the spectral function, J (oj) 

Efe \g k \ 2 S(uj-uj k ), 

- (1 — cosujt) doj. 
Similarly, the renormalized energy state is 
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duj. 



(12) 



(13) 



For an ohmic boson bath |43j , J (u>) = coo for u) < lu c . 
At high temperature with respect to its cutoff frequency, 
<j> (t) ~ r)t, where rj = aic/f3, and Eq = Eq + auj c . Gener- 
ally uj c is quite small compared to molecular energies, we 
thus ignore the energy shift, which is valid except when 
the noise strength is very large. This gives 



Gdna(E) = 



1 



E - E + 17 + ir/ 



(14) 



for the retarded Green's function. In this work, 77 — 
n/r n = 0,6.6 x 10- 3 ,6.6 x 10- 2 ,6.6 x 10-\6.6 eV for 
the timescales considered. For an interacting junction as 
given by the Hamiltonian in Eq. [8j the current is given 
by (using Eq. 4.114 in Ref. 



Iiji) 



2e 2 V 



7 



rn 



Et + { 1 + r,y m + fr+vY 



(15) 



The first and second terms represent precisely the first 
and second contribution in Eq. [3] respectively. However, 
as we have previously discussed, within this model calcu- 
lation, the liquid environment is allowed to form coherent 
interactions with the current-carrying electrons inside the 
junction. This results in the second term giving rise to 
orders of magnitude increase in the total current to values 
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that are unlikely in the present setting. In the junction, 
one has to consider also that the bosonic environment 
scatters the current-carrying electrons in all directions, 
including along the pore channel where they can be col- 
lected into the liquid. This effect is exacerbated by the 
fact that the environment both carries some longitudinal 
momentum and can act as a sink for electrons as well, due 
to the longitudinal bias. Therefore, on physical grounds, 
we assume that such processes occur which provide only 
the first contribution to the current in Eq. [3] Again, 
this is equivalent to assuming a two-probe noise model, 
as we have discussed previously. Under this assumption 
and for 7 <C Eq, the expression in Eq. 



15 becomes 



2e 2 V 7 2 



h 



El 



(16) 



i.e., the current for just a single structural distortion for 
linear response and weak coupling, and in the absence of 
inelastic processes that enhance the current. Note, that 
irrespective of this approximation, our main conclusions 
would be qualitatively unchanged. 

We know from above that the current acquires a distri- 
bution when structural distortions of the DNA are taken 
into account. Under the assumptions that went into 
Eq. [14] we can introduce these structural distortions by 
allowing Eq or 7 to acquire distributions. From Fig. [4] 
it is clear that the current distributions on a logarith- 
mic scale can be approximated as a Gaussian when no 
noise is present, which indicates that the coupling to the 
electrodes is controlling these distributions, as only the 
coupling fluctuates on an exponential scale. By assum- 
ing the coupling to both electrodes is identical, we miss 
structural distortions which bring the base into closer 
proximity to one electrode and farther from the other. 
However, this is unlikely to affect the essential physics. 

Now, let us calculate the distribution of 7's using the 
curve in Fig.|4]with no noise. Using the fact that the cur- 
rent distribution on a logarithmic scale is approximately 
a Gaussian, and that we are in a weak coupling regime 
(7 <C Eq), hx r y/'y m , where 7,„ is the maximum likelihood 
coupling strength, should also follow a Gaussian distri- 
bution. 



K ln 7/7m) 



1 



2na 2 



: exp 



(ln7/7 m ) 2 
2o* 



(17) 



with the standard deviation 07 = 07/2 w 0.45, where 
a 1 is the standard deviation of \nl/l m with 77 = and 
I m the maximum value. (Note that the relatively small 
standard deviation of the current distributions In I/I m , 
as seen m Figs, g] and [6j is a result of the control exerted 
by the transverse field [HI [20] . In the absence of such 
control the current distributions span several orders of 
magnitude and have considerable overlap |21j.) The max- 
imum, 7 m , appears at 6.8 x 10~ 4 eV, when Eq = leV, 
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FIG. 5: Current distributions of a model system for the Ade- 
nine nucleotide represented by a single energy level Eq. The 
current distribution on a logarithmic scale is taken to be 
Gaussian in similarity to Fig. [4] for no noise. As noise is 
turned on, at first the distribution does not change at all, but 
around 77 ~ Eo, where r\ = h/r n measures the strength of the 
noise, the distribution starts to shift. At larger r\, the peak 
of distribution shifts to lower values as r)~ 2 . The off-resonant 
tunneling, indicated by large Eq as measured from the Fermi 
level, "protects" the current distributions from noise. 



which is approximately the energy separation of Ade- 
nine's LUMO from gold's Fermi level pQ. We assume 
that the standard deviation of In 7/7™ does not change 
when we turn on the noise. The resulting current distri- 
butions are plotted in Fig. [5] 

Although we assume in our model that the distribu- 
tions stay Gaussian with the same standard deviation no 
matter what the noise strength, our model explains the 
key features found in our numerical simulations. The fact 
that the molecular energy levels are far away from the 
electrode Fermi level "protects" the distributions from 
this type of noise. This is represented by the term 
(Eq + rj 2 )^ 1 in the current (Eq. 16 1. The other features 



that appear, such as increased broadening and eventual 
multiple peak development, are not explained by our sim- 
ple model. These are due to multiple energy levels, Ei, of 
the fluctuating nucleotide junction, contributing to trans- 
port. The contribution from each reaches its turning 
point, rj « Ej,, at a different value of 77 and thus the 
single peak broadens and develops into multiple peaks. 

For the remainder of this paper we examine the role of 
transverse bias on the distributions for two different noise 
strengths (i.e., no noise and a timescale of r„ = 10~ 15 s). 
The results for the cases of 0.1 V and 1.0 V transverse 
biases are presented in Fig. [6j Previous work has shown 
that the transverse bias has a nonlinear effect on the 
mean of the distribution [50] . This is due to both a 
pulling effect of the backbone toward one electrode as 
the field is increased with consequent alignment of the 
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FIG. 6: Normalized current distributions for the four nu- 
cleotides at a transverse bias voltage of 1.0 V (top) and 0.1 
V (bottom). The solid lines correspond to an infinite noise 
timescale (no noise) and the dark dashed lines represent the 
distributions for t„ = 10~ 15 s, with the light dashed lines rep- 
resenting the bins used to produce the distributions. 



In conclusion, we have presented results combining 
molecular dynamics simulations with quantum mechan- 
ical current calculations including a model of noise gen- 
erated by the complex liquid environment in which the 
DNA translocation and interrogation takes place. We 
have shown that for reasonable timescales, e.g., down to 
1CP 15 s, the noise considered here will likely not affect the 
distinguishability of the current distributions obtained 
from measuring the transverse electronic current of the 
different DNA nucleotides. At extremely fast timescales, 
below 10 -15 s, the distributions are significantly altered, 
but this is beyond physically reasonable times for the ex- 
perimental system we are considering. We have also pro- 
posed a simple model system which provides insight into 
the physical mechanisms of noise and why the current dis- 
tributions are protected. This is due to the off-resonant 
nature of tunneling through the nucleotides and thus it 
is likely to be a general property of transport in organic 
molecules. While the distributions arc only mildly af- 
fected, we have shown that the type of noise we consider 
can potentially alter a single current count significantly, 
further supporting the notion that only a statistical study 
of the transverse currents can potentially distinguish the 
nucleotides. We finally note that while our study is done 
for a nanopore geometry, the results are applicable to 
other types of sequencing devices as well, such as the 
nanochannels of Refs. [2U HEj used in transverse elec- 
tronic measurements. 
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base toward the other electrode, and the steric effect of 
the alignment of the backbone with one of the electrodes. 
Therefore, while one can expect the mean current to be 
shifted to lower values with lower bias, the degree to 
which this occurs is not easy to determine a priori. This 
is especially true with the smaller base T. For this base, 
one cannot always expect perfect alignment at all times 
with the electrodes even in the presence of a stabilizing 
transverse field, further emphasizing the statistical na- 
ture of this problem. These effects can be seen in Fig. [6] 
In addition, one can see that all of the distributions are 
shifted slightly to lower current values due to noise, corre- 
sponding to an overall lowering of the current magnitude. 
However, the distributions themselves are very similar to 
the case of an infinite timescale (zero noise) . 
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The tight-binding Hamiltonian is constructed at 
every snapshot using the YAEHMOP package 
(http://yaehmop.sourceforge.net/), with the Fermi 
level also consistently calculated using this method. 
This condition entails that Ji + 12 = 0, where 1 and 2 are 
the two probes. This together with current conservation, 
It + h + Ii + Ja=0, yields It = —It, which is the current 
calculated in this paper. 
For simplicity, we set h — 1 in Eqs. 
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